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Figure 1. --Vapor pressure isotherms for temperatures of 72°F. and 
90 °F, show the moisture content m plotted as a function of the relative 
vapor pressure H. These curves are those which Spalt (9_) obtained for 
basswood for desorption conditions. 

Figure 2. --The parameter (3 is shown as a function of the moisture 
content m for temperatures of 72°?. and 90 °F. The ordinates of these 
curves are equal to the slopes of the curves in figure 1. 


Figure 3 . --Theoretical -drying curves for specimens with six different 
geometric forms: A, cylinder; E, infinite plane slab; C, infinite 

square rod; D, cube; E, sphere; F, semi-infinite solid. The ratio . 

is plotted as a function of the Fourier number 


OL 


Si. 


Figure 4. --The theoretical distribution of moisture across a plane 
slab at time t = ■ . 


Figure* 5 . --The relationship between the timeiag interval 'T and 

the specimen dimension r in which ' r /'7g is shown as a function of 

r/r on a log-log plot, qr was chosen as 1 hour for all specimens, 
o ° 

r was taken as 0.08&7 inch for the birch dowels and sawdust layers 
o 

and as 0.125 inch for the square basswood rods. 


Figure 6. --The quantity m - m for a ::-u;ure basswood rod 0.15 inch in 
diameter is plotted as a function of t/.e drying time for a temperature 
of 80°F. and a relative vapor pressure e 0.41. 



AN ANALYSIS OF TILE DRYING PROCESS IN FOREST FUEL MATERIAL 

George X. Byram 

Southern Forest Fire Laboratory, Southeastern Forest Experiment Station, 
U. S. Department of Agriculture, Macon, Georgia 

ABSTRACT 

It is assumed that the flow of moisture in forest fuels and 
other woody materials is determined by the gradient of a quantity g 
which is a function of some property, or properties, of the moisture 
content. There appears to be no preferred choice for this function, 
hence moisture transfer equations can be based on a number of equally 
valid definitions of g. The physical meaning .and dimensions of the 
mass conductivity will depend on the definition of g but the 

mass diffusivity <*: is independent of g. 

Simplified solutions to the transfer equations are expressed 
in terms of the ratio • The timeiag is a measure of the 

drying rate. It scales as the second power of the appropriate dimen- 
sion of the specimen, such as the half thickness, when the Biot num- 
% 

ber is large but scales as the first power of the characterizing 
dimension when the Biot number is small. Analytic solutions are 
not possible when properties are variable but the scaling relation- 
ships remain unchanged. 

The basic theory of the drying process throws considerable 
light on the complex interacting effects of air movement, radiation, 
and evaporation cooling on the crying rates of forest fuel material. 
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INTRODUCTION 


One of the most important factors influencing the intensity 
of forest fires is the amount of water in the burning fuel. Hence, 
any procedure for measuring or estimating the flammability of such 
fuels must give considerable weight to their moisture content and. 
its distribution. Forest fuels are a heterogeneous mixture of woody 
materials in which the individual components may vary greatly in 
size, shape, and arrangement. Intermingled with the nonliving fuel 
to a greater or lesser extent are green or living plants. They 
constitute another component of the. total forest fuel material. 

This paper is concerned primarily with the theoretical aspects 
of the loss of moisture in the nonliving fuels — especially in the 
moisture content range below the fiber saturation point . where the 
flammability of forest fuels increases rapidly with decreasing mois- 
ture content. The main' objective is to obtain a mathematical 
description of the drying process which can be readily compared 
with experimental data. 
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Among the numerous papers in the literature on drying, there 
are a considerable number which apply to the drying of wood or 
materials similar to wood. A few examples are the papers by Tuttle 
(11) , Sherwood (8) , Hougen, McCauley, and Marshall ( 5) , Stamm (10) , and 
Van Arsdel (12) . An extensive summary account of the drying of solids 
is given by Marshall and Friedman (7) . 

in lormuiat ing theories of the drying process, investigators 
have used different assumptions as to the nature of the gradient 
responsible for the flow of moisture in the drying material. One 
assumption is that a moisture gradient determines the flow, and some 
problems in diffusion are approached from this viewpoint. Another 
assumption is that the moisture flow is determined by the vapor pres- 
sure gradient. Van Arsdel (12) points out that, in the absence of a 
temperature gradient, the choice between these two assumptions is only 
a matter of convenience. However , other choices of gradient are 
equally valid and among the various possible gradients none can be 
designated as the preferred choice. This somewhat confusing point 

b 

can be clarified considerably by postulating that the moisture flow 
is determined by the gradient of a quantity g which '-is a function of 
some property, or combination of properties, of the moisture in the 
drying material. On the basis of this assumption one can then write 
the equation 
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2 &£l i s the 


wnere is tne rate of mass flow of moisture through an area A 

dt 

~J Q. 

perpendicular to the direction of flow, and id-2 is the gradient of 

dp 

g in the direction of flow. The quantity k is defined by this 

. & 

equation. It may be designated as the moisture or mass conductivity 
and corresponds to the thermal conductivity in heat transfer. Al- 
though k is numerically equal to the rate of moisture flow per unit 

S> 

of cross section area per unit of gradient of g, its physical mean- 
ing and dimensions will depend on the definition of g. 


MOISTURE TRANSFER EQUATIONS 

The fundamental equation of moisture flow based on equation 
(1) can be written as 



where C is the mass of water per unit volume, or the moisture con- 
centration in the drying material. This equation may be regarded 
as a* general form of Pick’s diffusion equation. The next step is to 
express g in terms of some property, cr\ properties, of the moisture 
content. One choice is the relationship g = P, which is equivalent 
to the assumption that the -moisture flow is determined by the vapor 
pressure gradient. ' If g = P, then equation (2) can be written as 




in which k has-been designated as k when g - P. 

O' Tj 

O JT 

At this point it is desirable to introduce 
quantity (3 which is defined by the equation 


(3 = 


j- 


the dimensionless 


where rn is moisture content and is defined as the mass of water per 

unit mass of fuel in the fuel -water mixture. The relative vapor 

pressure K is given by the relationship— ^ K = P/P where P is the 

s s 


ll H has values between 0 and 1 and m has values between 0 
and oo . For convenience they are usually multiplied by 100 an! 
expressed in percent j in which case Ii is known as relative humidity 
and m as percent moisture. 


saturation vapor pressure at temperature T. The defining equation 

ft 

for (3 shows that this quantity is the slope of a vapor pressure 
isotherm in which m is expressed as a function of K. 

The moisture concentration C may be expressed as 

C ~ C™ (4) 
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where |0 is the density^ of the fuel component of the fuel-water 


2/ The density is based on the bone dry fuel weight and 
the volume at moisture content m. 


'■■/Cure. If at any point in the drying substance the temperature T 
and density p are not changing with time, equation (4) can be 


differentiated and expressed successively by the series of equatione 

at 


St P \WWjj- S't & H Jj- e>"t~ 


_ p (3 ■§_£ 

~ p>~ S>Tr 

* 3 

Owing to shrinkage, the density p at a given point will decrease' 
with time as the drying progresses but the change is small and can 
be neglected. Equation (3) can now be written as 


7>fo$L &fa3§) + p6 2J? 


a 

s 




(5) 


( 6 ) 


This ‘equation describes the drying process even though and 
are variable* If k^ is constant, equation (6) takes the simplified 


form 


'al£ = 


(7) 
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where may be defined as the mass diffusivity and ie given by 
the relationship 


0 ( = 'hp ^ 

P(3 


( 8 ) 


One of the most effective choices for the function g is 
n 'zz-ym £ * 'If P is not changing with time then 

3 * 8 ; 


c? g __ D P^ 

^ ■£ " J 


(9) 


Combining equation (9) with the relationship := - and 


substituting the result for 




^•<L 


in the right member of equation (2) 


gives when written in terms of m 


j* 

^ A ' 


/° ^ 


( 10 ) 


If k is. constant and if there are no temperature gradients in the 
£ 

drying material^then equation (10) becomes 


y 

yyrj f p_j221 "' 2 ^ 2 . 

^ ^r z ^ 7 * 


^ L ~2svr? 

A.FZ ~ 

<7 5 


(ii) 
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A p 

<* = Plx-t 


( 12 ) 


Pairs of equations, analogous to equations (6) and (7), or (10) and 
(11), based on other choices of g will not be given. However, the 
relations between the k's and the & s for isothermal conditions are 

summarized for five different functions for g in table 1. All of 

> 2-1 
the cC 3 are identical and have the dimensions L t . Each of the 

k's is different, although for g * P and g « mP the two corres- 

6 

ponding k’s each have the dimensions of time. For simplicity the 

k for g *= mP is designated as k and all other k's are expressed 
s . g 

in terms of k^ . Table 1 illustrates the arbitrary nature of the 

gradient which is assumed to determine the flow of moisture. The 

procedure for finding the relationships between the k’s can be 

illustrated by determining the relationship between k and k . By 

g P 

letting g take the forms of g = nP and g » P it follows from 

S 

& 

equation (1) that for isothermal conditions 


2 sF ? ’ ^ 


Also for . isothermal conditions 


■A? ay? 
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Table 1 . — A summary of the k*s and <*.5 with their physical dimensions 
corresponding to five different choices for the function g 


. £ . 

k 

Dimensions 
of k 

2 -1 

cv (Dimensions are L t ) 

in? 

6 

j, 

t 

o( = 4 /* 

* 

? 


t 

* = £JL 

m 

Ar*. = £^5 

ml'±' 

u = 

C 

A = 

lA ~i~' 

* - -4# 

CP 

£ 

i = ^ 
cp P 

. aT'/ 3 * 

v = ^ 



i 
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Theref ore 


For a substance for vhich the vapor pressure isotherms are 

linear (Henry’s lav) (3 is- constant. In this case k vould be 

P 

directly proportional to k . If either k or k were constant, the 

-P 6 . 

other would also be constant. However, for most substances (3 is 
variable, so ordinarily both k's could not be constant. The vapor 
pressure isotherms shown in figure 1 for temperatures of 72°F. and 
90 °F. are those which Spalt (9_) obtained for desorption conditions 
for basswood. Figure 2 shows (3 plotted as a function of m for 
temperatures of 72°F. and 90°F. Ordinates of these curves were ob- 
tained by measuring the slopes of the curves in figure 1. Either 
curve shows that (3 can be. taken as approximately constant for 
values of m between 0.05 and 0.12 but increases rapidly a9 m approaches 
the fiber saturation. point . 

Figure 1. --Vapor pressure isotherms for temperatures of 
72°F. and 90°F. show the moisture content m plotted > 
as a function of the relative vapor pressure H . These 
curves are those which Spalt (9J) obtained for basswood 
for desorption conditions. 

Figure 2. --The parameter ( 3 is shown as a , function of 
the moisture content m for temperatures of 72°F,. and 
90°F. The ordinates of these curves are equal to the 
slopes of the curves in figure 1. 
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SOLUTIONS TO THE MOISTURE TRANSFER EQUATIONS 


When the fuel properties are variable, solutions to the 
moisture transfer equations have to be obtained by graphical or 
numerical methods. However, much can be learned about the general 
drying process from analytic solutions based on the assumption that 
the properties of the drying material, on certain combinations of 
properties, are constant. Such solutions are also useful approxi- 
mations even when the properties are not constant. 

A wide variety of solutions is available in the literature 
of heat and mass transfer. An important example is the solution for 
the infinite cylinder, which is significant for isolated individual 
materials like logs, dead limbs, and certain types of cured grass 
stems. Solutions for the infinite slab would apply to individual 
leaves, thin grass blades, isolated layers of bark and wood, and to 
a layer of duff- over a rocky surface impervious to the passage of 
moisture. Deep layers of duff and organic soils can best be treated 
as semi-infinite solids. For all of these cases the moisture trans- 
fer equations can be written in one dimensional form. 

Solutions based on equation (11) can be expressed directly in 

terms of the moisture content, so this equation is ordinarily the most 

suitable to use. Consider an infinite cylinder of woody material of 

radius r with a uniform initial moisture content of m which is 

o 


J 
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drying in an atmosphere of constant temperature and constant relative 
vapor pressure for which the equilibrium moisture content is m^. If 
m is the moisture content at any radial distance x from the cylinder 
axis at time t, equation (11) in cylindrical coordinates becomes 


JL - -i- 2^1 


The initial conditions are 


m m ** m - m 
e o e 


when t ■ 0 


The boundary conditions -are 






at x ■ 0 


and 


* P I f \ 


7 ^ 7 - 


at x « r 


where h is the surface transfer coefficient* 
g 



3 / 

The solution to equation (13) is— 


cO 

~ k-A %' 

sn - 1 


3/ Details of solutions for objects of various shapes are 
•given in works on diffusion and heat transfer such as those of Crank 
(3_) , Gied.t (5_) , Eckert and Drake ( 40 , and Carslaw and Jaeger (2), 

As in heat transfer theory, the dimensionless groups /C^l 7 ’ and 

h r/k may be designated as the Fourier number and Biot number, 

s s 

respectively. 



where ** %^r and ^ n ’ represents roots of the eigenfunction equa- 
tion 

(a^j = (X5) 

J is«a Bessel function of the first kind and zero order and J, is a 
o 1 

Bessel function of the first kind and first order. 

In experimental work, the average moisture content is much 
easier to determine than the moisture content at a given point with- 
in a specimen, so it is desirable to have solutions to the equations 
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expressed in terms of the average moisture. If ra is the average 
moisture content throughout the cylinder then 


s?71 — 




V 

Substituting in this equation the value of m - m^ given by equation (14) 
and integrating gives 


n y - 


OP *!■> y 

sn c t 


-al 1 


(16) 


The value of in equations (14) and (16) will depend on the 

ratio h r/k in equation (15). The surface transfer coefficient h 
g g g 

is defined by the general equation 




(17) 


where is the rate of mass flow of moisture through a surface of 

area ‘A, is the value of g at the surface of the specimen^ and is 
the equilibrium value of g for an environmental relative vapor pressure 


H^. For g *■ m? e equation (17) becomes 


^r= 


(18) 
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Equation (16) takes on a much simpler form if the resistance 

to the flow of moisture through the surface 1 b negligible compared 

to the internal resistance to moisture flow. In this case the Biot 

number h^r/k^ is very large. In equations (14) and (16), J^N^) — 

as h r/k * oo . Equation (16) then takes the limiting form 

g g 


Tn* 



/77-i 


-AC fa*/''* 3 ) 

e 


o 



in which N, , FL ........ N are the successive zeros of J (N ) . 

1 ' * ■ 2 * n o n 

Equation (19) is represented by curve A in figure 3, in which « / ^ 71 ~ 
is shown on a semi-log chart as a function of the Fourier number 


Figure 3. --Theoretical drying curves for specimens with 
six different geometric forms: A, cylinder; B, infinite 

plane slab; C, infinite square rod; D, cube; E, sphere; 

F, semi-infinite solid. The ratio i.s 

plotted as a function of the Fourier number- . 



With the exception of the brief initial period vhen the moisture 
distribution is becoming established, this curve is linear. This 
is a result of the rapid convergence of the series in equation (19). 
For values of greater than 0.10, only the first term need be 

retained^ and equation (19) can be written approximately as 

__ .. -a 

nm — £ 

. w 


( 20 ) 


For an infinite plane slab, equation (11) reduces to 

^ "V/? __ ji 

- PL V *tr 


in which x is measured from the slab center. The initial conditions 


and boundary conditions are the same as for equation (13). For a 


3/ 


slab of thickness 2r, the solution is— 


oo 


s77j- syrjjt 




— A. 




✓7 */ 


(\/ *“f- -f /V^ 


-N?(*+A*) 

1 Cjr* (a '~r9f) 


( 21 ) 


where * and represents the roots given by the 


n 

equation 


= 4 ^ = 

7 3 * 


( 22 ) 
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When expressed in terms of m - , the solution for the 


oO 


rh 7 —rrru. 

“'^7 J e 








sn 


( 23 ) 


When h r/k is very large, equation (23) also takes a simpler 
g S 

form. As h r/k — 

8 g 

Equation (23) then becomes 


;> OO , cotK 0, and N . 

^ n * n 


00 


/ 


S Yri- 'TV* x _jl_ £ 




•^ur^u 


Hi, 


± 

== TT' 


I w ' — i p ^ / j«. - r U— « • • • 

Je e T* 3-*^ -r 


(24) 


which is represented by curve B in figure 3* As in the cylindrical 
case, the series in equation (24) also converges very rapidly and 
the straight.line portion of curve B represents only the first term, 

k 

Equation . (24) then reduces to 


i Z>'V, 


^v>7 — syyjjt 

syyij-sryi* 


TT' 


(25) 
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3 / 

By the use of the product method*- , solutions in one 
dimension can be applied to objects of numerous shapes. An ex- 
ample is the infinitely long square rod. In this case it can be 
shown that the solution is the square of the right member of equation 
(23) or equation (24) for the infinite plane slab. If only one term 
in the resulting series iB retained, the equation analogous to equa- 
tion (25) 1 b ' ' > - 



where 2r is the diameter of the rod. Similarly, the equation for a 
cube of diameter 2r is 


rm - 



-3 (?) 


( 27 ) 


In figure 3^ equations (26) and (27) are represented by the straight 
line portions of curves C and D, respectively. 

Curve E in figure 3 is the drying curve for a sphere of radius 
r. The straight line portion of the curve represents the equation 





-T 

e 
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Equation (20) as veil as equations (25) through (28) can 


all be expressed by the single equation 


s yrl 




( 29 ) 


in which ie a quantity which has the dimensions of time. It 

represents the combination of factors which comprise the coeffi- 
cients of t in the different drying equations. Physically qr is 
equal to the time required for the value of fn - at some arbi- 
trarily chosen point on the linear portion of any of the curves in 
figure 3 to drop to 1/e of the value at the beginning of the interval. 
It is defined as the timelag and is a very useful parameter for 
comparing the drying characteristics of different forest fuel 
materials. The dimensionless constant K depends on the shape of 
the drying specimen. Geometrically it represents the intersection 
of the extension of the straight line portion of the curves in 
figure 3 with the t * 0 ordinate. For example, in curve B, the 
value of K is <5/V- , . Equation (29) can also be written as 
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in which m - m ■* K(m - m ) and represents an adjusted or 
o e o e r 


equivalent value of m - for t *= 0. 

An additional case of considerable practical significance 
is the semi-infinite solid which is approximated by a deep layer 
of duff or organic soil. Throughout a surface layer in which the 
depth r is small compared to the total depth, m - can.be expressed 

as an integral of the error function. If constant properties are 
assumed and If Is very large, the equation for m - m^ can be 


written as 




/-W7 — — — 

~ /x T'7rf 



_ u 




— u. 


'o 'o 

Placing 6? in a power series form and integrating term by 

term gives 


✓^77 — SVK . 


£—\i - ..... 7 

Xfn [' IrTT. ^ j-jT.Zi '/ ■ 7 ■ J 


(31) 


where This equation is represented by curve F 

1 ' 

in figure 3. In this case the logarithm of To - does not become 
linear with increasing values of as in all of the other 

cases but decreases at a diminishing rate. This is caused by mois- 
ture traveling upward from deeper layers in the semi-infinite material. 
The thickness of the surface layer is taken as r rather than 2r be- 
cause, unlike the slab, it can lose moisture from only one side. 
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DRYING RATES 


The Effect of Temperature and Relative Vapor Pressure 


To examine the effect of temperature and relative vapor 
pressure on the drying rate, it is somewhat simpler to work with 
vapor pressure gradients rather than with gradients of mP . Tak- 
ing the infinite slab of thickness 2r as an example, the solution 
of equation (7) in one dimension is 


P~ 2\ 3 ' 

/§“(? / ' /V Tf ‘f-O'v/t, 

u m = / 





where P is the vapor pressure in the drying material at a distance x 
from the slab center, P^ the yapor pressure at infinity, and P q the 
vapor pressure within the specimen at time t ** 0 when the moisture 
content m was uniform. 

The form of the right members of equations (21) and (32) are 

identical but N^is not the same in both cases. In equation (21) 

N is a function of h r/k and in equation (32) it is a function of 
n g g 

^ • The surface transfer coefficient h p is defined 

by. equation (17) with P i replacing ^ , P^ replacing 

and replacing . 

The outward rate of mass flow of moisture across the surface 


is 





where 


d K 

dt 


now represents the mass flow per unit of surface area 


and ^ va P or Pressure gradient at the surface. 

Differentiating equation. (32) and letting x » r gives 



4^ P 6 /^) 


P-P 

R-e ~~ 3 


SSL 



ti 


^77 = / 


( 34 ) 


(35) 


where is the vapor pressure at the slab center. Combining 
equations (33), (34), and (35) gives for the outward moisture flux 

fjf- “ # ft - % 


(36) 


where is the sum of the series in equation (34) and the sum 
of the series in equation (35). A drying rate ratio R may be de- 
fined as 


R 


diM , ,dH ' 
dt ' Mt j 
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dM 

in which — is the outward flux of moisture for vapor pressures 

. dt , 

of P and P^ and the ^ or vapor pressure P^* and P^ 1 « 

The drying rate ratio makes it possible to compare the drying rate 
of two identical specimens in environments with different tempera- 
tures and different relative vapor pressures. 

Using equation (36), the drying rate ratio can be written 


as 



1 



If is chosen so that this group is the same for both speci- 
mens, then ■ S^ 1 and . R can then be written in the 

alternate forms 

r- !kJz - tLztk 
. /?'-&' ni~ & ~~R' 

where P and P 1 are the saturation vapor pressures for temperatures 

S 6 

T and T 1 in the two environments, H^and are the corresponding 
relative vapor pressures in the environments, and and IT 1 are 
the relative vapor pressures at the slab centers. 

If the moisture content at the center of the two slabs is • 
high,) then P^ and P^’ are nearly equal to the saturation vapor 
pressures P and P 1 , respectively. Equation (37) can then be 

S 6 



written in the approximate form 


o (7 / tj ) B 

- (Pfr 


( 38 ) 


which can be used for estimating R when the central moisture content 
m is 0.20 or more and if H and H’ are 0.40 or less. For smaller 

C Jc ■ JC 

values of m and higher values of H and Hi, equation (37) should be 
. c ■* -* 

used. The values of H and H 1 corresponding to m and m * can be 

c c c c 

read from the vapor pressure isotherms. The central moisture can be 
determined from the average moisture which for the slab is given by 
the . equation 

rm- ^]a = -Sf ~ y ^ r U) 


provided that h r/k is large and the drying periods are long enough 
S S ' 

for the moisture distribution to become established. This distribu- 
tion begins to approach its final form at relatively small values of 

as shown in figure 4. This diagram is a plot of equation 


Figure 4. --The theoretical distribution of moisture 
across a plane slab at time t * */i.. 


(21) for t ” ^ /^L i or a ~ ^ f TT^' > and h r/k ” c 0 

and represents the theoretical distribution of 2 21112^ across the 

syrijr ^ 

slab. 
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It can be shown that equations ( 37 ) and (38) apply to speci- 
mens of other shapes such as cubes, spheres, and rods* However, for 
a finite surface layer of. a semi-infinite material, the drying rate 
ratio becomes 




'/JL 


For two identical specimens drying in environments for which 
the specimen temperatures are T and T 1 the ratio of their timelags 
will be 


Vr 1 = Vr 


This ratio should not be confused with the drying rate ratio R in 
equation (37) although both are a measure of relative drying- rates. 

R represents the ratio of the rates when either rate is expressed 
in units of mass, of water per unit time; ^ * re P re6ents the ratio 

of the reciprocal of the rates of change of which, with 

the exception of the semi-infinite case, are constant when the dry- 
ing curves (figure 3) become linear. The magnitude of t/ ' does 

k * » 

not depend on the relative vapor pressure-. 


The Effect of Fuel Size or Thickness of Fuel Layer on the Drying Rate 

Equations such as (14) and (21) express m - as a function 
of three dimensionless groups , and x/r. ' In equations 




for m - m , such as (16) and (23), only the groups and cc£ /si?~ 
appear because x/r *= 1. It is not readily apparent from any of 
these equations Just how the time required. for m - m^,or ® ” m e > 
to decrease to some given fraction of m ■- depends on the dimen- 
sion r because the values of also depend on r. However, , 

K , . . . . . N approach constant values in the limiting cases 
2 n 

for which h r/k > co and h r/k >• 0, Examples of the 

g g g g 

first case are represented by equations (19) and (24) ? in which only 
the group appears as a variable in the right member of the 

equations. In this case it is obvious that the time required for 

_ 2 
m - m to drop to a given fraction of m - m Varies as r . The 
e r o e 

same relationship holds for specimens of other shapes, including 
a finite surface layer of thickness r which is part of a semi- 
infinite material. The second power relationship also applies to 
m - at any given value of x/r in equations (14) and (21) as 

h r /k — ^ co except at the surface where x/r ** 1 and 

g g ' 

m - m^ » 0, • The same results are obtained with equations based on * 

the other functions for g given in tab lev 1. 

The form which the equations for in - m .takes when h r/k — > 0 

eg g^" 

can be most easily determined from equation (23). — > fa - 0 77 

h r/k — > 0 and all terms in the series vanish except the first for 
g g J 
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■which n =1. This term is the product of the term 0 
and the indeterminate fraction 






L 


/V, fort 


which approaches 1/2 as N 


0. From equation (22) 


/V N, = h# 


Since tan N 


1 


N x as N 1 « 


9 


0, it follows that N 


h r/k for 
g g 


very small values of h r/k • In this case equation (23) thus takes 

8 g 

the . limiting form 


srri-srru _ p J V 

srnji. c 


=r e 


—4 ^ 4 //°^ 


The time required for m « m^ to reach a given fraction of 

thus varies as r when h r/k is very small. Equation (21) also 

g g 

approaches equation (39) as h r/k ^*0. This means m = m for all 

g g 

values of r and moisture gradients in the drying material vanish in 
this limiting case. 

From its definition the timelag qr for specimens of different 


shapes can be expressed in the forms 

■7 


* SO 



For large values of h r/k . N. is a constant whose value depends 

g g 1 

on the specimen shape. ; For very small values of h r/k , N. becomes 

£ g 1 



3 


^3 


where B is a dimensionless constant which also depends on the speci- 
men shape. When h r/k is very small, the equation for 7" then 

S S 

becomes 

(40) 

Thus there is a second power relationship between qr and r when 
the Biot number is large and a first power relationship when the 
Biot number is small. 

The derivation of the relationships between qr and r has 

been based on the analytic solutions to the moisture transfer equa- 

tions. In turn, such solutions, were possible only because of certain 

* 0 

simplifying assumptions including that, of constant k . The question 

S 

then arises as to whether the relationships between 7- and r hold 
when k^ is variable. It can be shown that they do by means of a 
similarity argument. 

, * 

In equation (10), let x, y, z, and t be replaced with dimen- 
sionless X, Y, Z, and t’ such that 

X * x/r, Y * y/r, Z = z/r, and t ‘ - tjqr 



- 28 - 



where r is some appropriate dimension of the specimen such as its 
radius or half thickness and T is the tiraelag for the specimen 
with the dimension r. Also let be variable and given by the 


equation 




where f(m) is a dimensionless function of m and k 1 is constant 

8 

and represents the value of for some given Value of m. Equa- 
tion (10) can then be written as 

b (fly?) c>(+1sy) 9^ 1km - km, 

kFx~ ^ T ^ 


(41) 


where cK mm 's J jp • The dimensionless group o( *T /ki ^ 
corresponds to the Fourier number 0( "L /. A which appeared in the 
analytic solutions* Corresponding to the Biot number is the group 
/y Ji f which appears when the . surface boundary condition is 

$ 3 


9 

written in dimensionless form. However, as h r/k 1 
. g g 


QQ 


the solution of equation (41) in terras of .’m -‘m will be a function 
only of (<X / 'r/>S)~t t • For specimens of the same shape, the 


group 


o< 1 T /x- 


is constant. Hence^for a given value of o( 


T 


must vary as r 
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The relationship between T and r when h r/k — — v> 0 can 

g g 

be determined in the general case from equation (18). k may be 

g 

variable within the specimen, but since moisture gradients approach 

zero as h r/k ■ — ^ 0. the moisture content m can be considered con- 
g g 

stant throughout the specimen or m ® m. These conditions would be 
approximated by a block of wood coated with a paint film with a 
high resistance to the passage of moisture. If A is the total sur- 
face area of the specimen and V its volume, equation (18) can be 
written as . 

(oVj £e? = -£/4 A 

The ratio V/A can be written as Br where B is the same dimensionless 
shape factor that appears in equation (40). Equation (42) can then 
be written as 

d (m- )_- !?K _ _ J__ 

which is identical to equation (40). 



Preliminary Drying Rate Measurements 

Although the details of the experimental work will not be 
given in this paper, some of the preliminary results of drying 
experiments will be briefly summarized to show the relationship 
between '7* and the specimen dimension r for several different 
materials. Figure 5 shows a dimensionless log-log plot of rr /% 
and r/r^ in which r f is the timelag for a specimen temperature 

Figure 5. --The relationship between the timelag interval 

T and the specimen dimension r in which /rj~ is 

shown as a function of r/r^ on a log-log plot, 'jr was 

chosen as 1 hour for all specimens, r was taken as 

o 

0.0867 inch for the birch dowels and sawdust layers and 
as 0.125 inch for, the square basswood rods. 

of 80°F. and r is the radius or half thickness depending on the 
type of specimen. The values of *T C and r^ were chosen so that the 
average position of the data on the chart would be approximately 
the same for the three different kinds of specimens. These were 
square basswood rods which ranged from 0.07 inch to 3.70 inches in 
diameter, round birch rods from 0.063 inch to 1.00 inch in dia- 
meter, and loosely packed sawdust layers 0.27 inch to 4.00 
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inches thick. The sawdust layers were in metal trays and could 
lose moisture from only one sidejSo the effective range of thick- 
ness was double the preceding range of values. The line through the 
data represents the theoretical second power scaling relationship 

between 7 - and r which should exist when h r/k is large. With the 

S £ 

exception of the points representing the smallest birch rod and the 
smallest basswood rod (and also the points for the two thinnest saw- 
dust layers), the agreement of the data with the theoretical line is 
good. Possibly for the smaller specimens^ h^r/k^ is small enough to 
cause deviation from the second power relationship between 7 ' and r. 
It may be that solid woody materials with r less than 0.02 or 0.03- 


inches, or sawdust layers with r less than 0.1 inches, are in the 
region where 7 " varies as the first power of r. However, this has 
.not yet been determined experimentally. 

Figure 6 is an example of the drying curves^ each of which yield 
one value for 7 * .In the basic curve, the quantity m - m^ is 


Figure 6 . --The quantity m - ra^ for a square basswood rod 
0.15 inch in diameter is plotted as a function of the 
drying time for a temperature of 80°F. and .a relative 
vapor pressure of 0.41. 



plotted as a function. of the drying time t on a semi-log chart. 

With the exception of the first few measurements, the points tend 

to fall on a straight line, the slope of which determines the time- 

lag 7" . The intersection of this line with the t ■ 0 ordinate 

gives the value of ^which appears in equation (30) . 

The value of k for any given specimen can be estimated 
£ 

from the equation 

~ AT,* Tp s 

Excluding the smallest specimen in each group, the mean 

values of k for' the birch and basswood rods was 2.34 x 10 ^ 

8 

seconds and 3.12 x 10.^ seconds, respectively. For the sawdust 

layers the mean value of k was much greater, about 17.4 x 10 ^ 

& 

seconds with the thinnest layer excluded. Apparently vapor diffusion 

in the air spaces of the sawdust made it a much better conductor of 

moisture than solid wood substance. This effect should be present 
« 

even in a greater degree in sane types of forest fuels such as 
layers or beds of conifer needles. 

The data for figure 6 were obtained by weighing the basswood 
specimens (six square rods 12 inches long and 0.15 inch in diameter) 
at about 10-minute intervals while they dried in a room of nearly 
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constant temperature and constant relative vapor pressure. At the 
end of the drying period^ the rods were ovendried and the m - 
values computed. The value of m^ was taken as 0,09 throughout the 
period. 

Possible Effects of Air Movement, Radiation, and Evaporation Coolin: 


There are complex interactions between air movement, solar 

radiation, and evaporation cooling in their combined effect on the 

rate of drying of forest fuels. With the exception of the early 

stages of the drying period, evaporation cooling should have little 

effect on the drying rate if h r/k is large. In the absence of 

g g 

radiation, wind should also have but little effect on the drying 

/ 

rate of solid material if h r/k is large. Byram ( l-2~) showed that 

S S 

specimens in bright sunlight will have a lower drying rate if there 

is wind than if there is no wind. In this case, the convective 

cooling of an airstream partially offsets the temperature rise due 

to radiation and results in a smaller value of P • If h r/k was 

s g g 

very small, the combined effect of wind, sun, and evaporation 

would be determined by the product h P in equation (43). Air 

g 8 

movement would increase h and . radiation would increase P . On 

g 8 

the other hand evaporation cooling and convective cooling (in the 


presence of intense radiation) would tend to limit the increase of 
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Although the relationships have not been determined experi- 
mentally, air movement should have a greater effect on the drying 
of porous materials 6uch as beds of conifer needles or hardwood 
leaves, than for solid material. The transport of water vapor out 
of the interior spaces by penetrating air currents should result 

in a marked increase in k . 

S 


CONCLUSION 

Depending on the nature of the moisture flow problem, some 
choices of the gradient, or "driving force," which induces the 
moisture flow, are more suitable than others. However, in a mathe- 
matical sense there appears to be no preferred choice. 

Analytic solutions to the moisture transfer equations, based 
on the assumption of constant properties, give approximate but useful 
descriptions of the drying process for a variety of forest fuel 
materials. In addition they result in a better insight to some com- 
plex features of moisture flow, such as -those involving surface 
transfer phenomena. However, some of the most significant relation- 
ships, such as the first and second power scaling relationships be- 
tween 7* and r, can also be deduced by other methods. These do not 
require the assumption of constant properties, nor do they require 
a solution of the differential equations. Hither approach leads 


- 35 - 



to basic relationships which are helpful in understanding and 
predicting the interacting effects of wind, solar radiation, and 
evaporation cooling in the drying of forest fuels. 
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